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ABSTRACT 



Accretion of minor satellites has been postulated as the most likely mechanism to explain 
the significant size evolution of the massive galaxies over cosmic time. Using a sample of 
629 massive (M stal ~ lO n M ) galaxies from the near-infrared Palomar/DEEP-2 survey, we 
explore which fraction of these objects has satellites with 0.01 < M sat /M centra i < 1 (1:100) 
up to z — 1 and which fraction has satellites with 0.1 < M sat /M centra i < 1 (1:10) up to z — 2 
within a projected radial distance of 100 kpc. We find that the fraction of massive galaxies 
with satellites, after the background correction, remains basically constant and close to 30% 
for satellites with a mass ratio down to 1:100 up to z = 1, and ~15% for satellites with a 
1:10 mass ratio up to z — 2. The family of spheroid-like massive galaxies presents a 2-3 times 
larger fraction of objects with satellites than the group of disk-like massive galaxies. A crude 
estimation of the number of 1 :3 mergers a massive spheroid-like galaxy experiences since z~2 
is around 2. For a disk-like galaxy this number decreases to ~1. 
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The relevance of major mergers 
the size increase of the massive 

the l a st -11 Gyr (e g.. iDaddi et alj 120051: iTrujillo et al.l 120061 
120071 : lLonghetti etal] 120071 : iBuitrago et afl l2008h has been dis- 
favor e d observationally (e.g.. | Bundv et alj 120091 : Ide Ravel et al.l 
2009: lLopez -Saniuan et al. I l201(iir This has left room to a grow- 



ing consensus that the strong size evolution observed among the 
massive galaxies is mainly dominated by the continuous accretion 
of minor satellites. However, all the observational evidences com- 
piled so far suggesting that the minor merging is the main route of 
galaxy size growth it is only indirect. The observations that favor 
the minor merging scenario are: a) a progressive build-up of the 
envel o pes of the massiv e galax ie s with cosmic time jHopkins et alj 
20091: iBezanson et alj 120091 ; Ivan Dokkum & BrammeJ l20ld : 
Carrasco et alj l2010h and b) a mild decrease in the veloc- 



ity dispersion of these galaxies (e.g., ICenarro & Truiillol 120091: 
Cappe llari et alJl2009l:lMartinez-Manso et al.ll201ll : lNewman et all 
2010; van de Sande et al.l201 lh . Both phenomena agree with a pro- 



cess that do not affect dramatically their inner regions. Recently, an 
extra evidence supporting the merging scenario has been stated: 
the size evolution of the massive galaxie s is not linked to th e age 
of the stellar population of the galaxies (Truiill oet alj|201lh . All 
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these observations disf avor the puffing-up mechanism proposed by 
iFan et al. I2008. l20lti) , where galaxies grow by the expulsion of 
gas by the AGN activity, and give support to the minor merging 
hypothesis. 

On the theoretical side, N-body c osmological simulations as 
well as semiana l ytical models (e .g., iKhochfar & Burkertl |2006l : 
iNaab et alj|2009l ; [Oser et alj|201 ll) show that the expected accre- 
tion rate of satellites should be able to produce a significant in- 
crease in the size of the galaxies while at the same time chang- 
ing the v elocity dispersion only mi ldly. Estimates of the merger 
rate (e.g.. Lopez-Sani uan et alJuOllI) using observations are, how- 
ever, not straightforward due to the large uncertainties in the de- 
termination of the merging timescales. Nevertheless, a more di- 
rect way of confronting simulations with observations and, con- 
sequently, probing the minor merging scenario is to measure the 
frequency of satellites found around massive galaxi es and quantify 
how this fraction changes with cosmic time (e.g.. iNewman et al] 
1201 lh . Several pa pers have calculated this nu mber in the nearby 
Universe (see e.g., IChenll2008l ; rLiu et alj|201 ]]). These works show 
that ~12% of the massive galaxies have at least a satellite with a 
stellar mass 0. 1 < M sat /M centra i < 1 within a projected radius of 
100 kpc. These numbers are in very nice a greement with expec- 
tation s from ACDM simula tion s (see e.g.. iBov lan-Kolchi n et alj 
l20ld) . lde Ravel et alj j20TTh and lNierenberg et alj J201 lh have ex- 
plored the evolution of the fraction of galaxies with satellites up 
to z ~ 1 but using mostly samples of central galaxies less mas- 



©2011 RAS 



2 E. Mdrmol-Queralto et al. 



sive than 10 u M o . In this paper we concentrate on the most mas- 
sive galaxies and we expand on the previous analysis exploring 
the fraction of galaxies with satellites up to z ~ 2. To reach our 
goal we use a large and complete sample of massive galaxies up to 
z = 2 from iTruiillo et al, I d2007h . We probe two different redshift 
ranges: up to z = 1, we explore the fraction of massive galaxies 
with satellites having 0.01 < M sat /M ccnlra i < 1, and up to z = 2, the 
fraction of massive galaxies with satellites within the mass range 
0.1 < M sat /M centra , < 1. 

This letter is structured as follows. In Section [4] we describe 
our sample of massive galaxies and the photometric catalog we 
have used for identifying their satellites. Our criteria for select- 
ing satellites as well as our background estimation methods are 
explained in Section [3] Finally, our results are presented in Sec- 
tion|4] and a discussion of our findings in provided in Section[5] In 
this paper we adopt a standard ACDM cosmology, with fl m = 0.3, 
D A = 0.7 and H = 70 km s~' Mpc~'. 



2 THE DATA 

To analyse the evolution with redshift of the fraction of massive 
galaxies having satellites, we have used as the reference catalog 
for the c entral galaxie s the c ompilation of massive objects pub- 
lished in ITruiillo etall d2007h (hereafter T07). This is a homoge- 
neous and large collection of massive galaxies since z = 2. Briefly, 
the sample consists on a total of 831 massive (M star >10"M o ) 
galaxies (of which 35 where identified as AGN and not used 
onwards) over 710 arcmin 2 in the Extended Groth Strip (EGS). 
These objects were K,-band selected in the Pal omar Observatory 
Wide-Field Infrared ( POWIR)/DEEP-2 survey dBundv et alj|2006l : 
IConselice et al.l |2007[). In total. 372 galaxies have spectroscopic 
redshifts dDavis et al.l l2003). whereas the remaining redshifts were 
obtained photometrically using B, R and I bands from the CFHT 
3. 6m- telescope, F606W and F814W from the Hubble Space Tele- 
scope and J and K, from the Palomar 5-m telescope. Stellar masses 
and other d erived photome tric parameters were estimated using 
a Chabrier (Chabrier 2003) initial mass function (IMF). T07 es- 
timated (circ ularized) half-light radius (r e ) and Sersic indices n 
dSersid 19681) for all the galaxies in our sample. 

To compile the sample of the satellite galaxies around our 
massive objects we have used the EGS IRAC-selected galaxy 
sample from the R ainbow Cos mological Databasj^] publis hed by 
iBarro et alj d2011al) (see also IPerez-Gonzalez et al.l 12008). This 
database covers an area of 1728 arcmin 2 centred on the EGS and 
provides spectral energy distributions (SEDs) ranging from the UV 
to the MIR regime plus wel l-calibrated and rel iable photometric 
redshifts and stellar masses dBarro et a l. 2011b). Around 10% of 
the galaxies in the Rainbow catalogue have spectroscopic redshifts. 
From the Rainbow database we have selected all the galaxies with 
z < 2.2 and an estimated stellar mass 10 s M < M < 10 12 M o . A to- 
tal of ~55000 objects were selected in the EGS area following these 
criteria. We refer to this resulting sample as the Rainbow catalog. 

The sample of massive galaxies as well as the Rainbow sample 
were cross-correlated using a 1.0" search radius to create a sam- 
ple of central galaxies identified in both catalogs. All the massive 
galaxies in T07 were found in the Rainbow database. The aver- 
age difference between the photometric redshifts for the massive 
galaxies in both samples is ~10%. The average stellar mass of our 
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Figure 1. Stellar mass vs. redshift for the massive galaxies analised in this 
work. Galaxies with spectroscopic redshifts are plotted in red, while galax- 
ies with photometric estimates are plotted in blue. The left panel shows the 
distribution of the massive galaxies used selected for exploring the frac- 
tion of galaxies with 1:10 satellites up to z = 2. The right panel shows 
the massive galaxies used in the study of 1:100 satellites up to z = 1. The 
solid green lines illustrate the stellar mass 75% completeness limit o f the 
Rainbow database for the redshift ranges given in Barr o et al] j201 lbl) . The 
magenta dashed lines the stellar mass cut used in this paper for the different 
subsamples. 



massive sample according to the Rainbow dataset is 0.9x10' 'M Q 
and 1.7xl0"M o according to T07. 

In order to build a sample of central galaxies with the best esti- 
mations of redshifts and stellar masses we have applied the follow- 
ing rules: a) if a central galaxy has a spectroscopic redshift determi- 
nation in Rainbow (348 objects), we have used this redshift plus the 
stellar mass inferred in that catalog for these two quantities. Among 
these galaxies, there were 8 objects with spectroscopic redshifts in 
both samples with high discrepancies in the stellar mass estimations 
from both catalogs. We reject from our sample such dubious cases, 
b) If no spectroscopic redshift is found on the Rainbow database 
but on the T07's sample (37 galaxies) we use the values of redshift 
and stellar masses from that catalog, c) Finally, if no spectroscopic 
redshift is found in any of the two catalogs we have used only those 
objects where the photometric redshift determination is robust (317 
objects). This means that we have co mpared the two inde pendent 
photo-z estimations found in T07 and lBarro et al] d201 lab and we 
have only taken those objects where the photometric redshifts dis- 
agree less than Az pho t = 0.070 for 0.0 < z < 0.5, Az pnol = 0.061 
for 0.5 < z < 1.0, and Az phol = 0.083 for 1.0 < z < 2.5 (typical 
quality of the photo-z's in the Ra inbow catalog in EGS obtained by 
comparing them with spec-z's. IBarro et al]|201 lbh . This removes 
94 galaxies. For consistency with the sample of satellite galaxies, 
for these 317 objects we take the stellar masses and photometric 
redshift from the Rainbow catalog. After this selection, the num- 
ber of objects in the final sample of massive galaxies is 694, of 
which 317 have photometric redshifts from Rainbow and 377 have 
spectroscopic redshifts (340 from the Rainbow catalog and 37 from 
T07). 

A final cut in the number of galaxies of our main sample is 
required to assure that the fraction of galaxies with satellites along 
our explored redshift range is not biased by the stellar mass com- 
pleteness limit of the Rainbow database. The stellar mass limit 
(75% co mplete) of the Rainbow data base at each redshift is pro- 
vided in Perez-Gonz alez et all d2008L see their Fig. 4). In the red- 
shift range < z < 2 we have selected only those massive galaxies 
whose stellar masses are 10 times larger than the completeness limit 
at each redshift. There are 629 galaxies (with a mean stellar mass 
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Figure 2. : Left panel: ACS color image of the massive galaxy IRAC123359 (in the centre) at z = 1.17 with a satellite galaxy (IRAC123191-1) that meets the 
selection criteria used in this work enclosed by a white circle. A circle of radius 100 kpc is plotted with a dashed line. Central and right panels: Spectral energy 
distributions for both the massive (central panel) and the satellite (right) galaxies. These panels also include the redshift and stellar mass estimates according 
to the Rainbow database. 



of M = 1.3 x 10 U M G for this sample) that meet these criteria. On 
doing that we are secure that we can explore within the Rainbow 
catalog satellites down to 1:10 mass ratio of the central galaxy in 
the range < z < 2. For the same reason, this exercise is done 
up to z = 1 but this time only selecting those central galaxies with 
a stellar mass 100 times above the mass limit. This cut leaves us 
in the redshift range < z < 1 with 194 massive galaxies (with a 
mean stellar mass of M = 1.7 x 10" M Q . The stellar masses and the 
redshifts for the central galaxies studied in this work are illustrated 
in Fig.Q] 



3 SELECTION CRITERIA 

To identify the satellite galaxies around our central objects we have 
applied the following procedure: (1) we identify all the galaxies in 
the Rainbow catalog which are within a projected radial distance 
to our central galaxies of -K scil rch=100 kpc (corresponding to 0.3 and 
0.2 arcmin for z = 0.5 and z = 2, respectively); (2) the difference 
between their photometric redshifts and the redshift of the central 
galaxies is lower than the lcr uncertainty in the estimation of the 
photometric redshifts of the Rainbow database (i.e., Az p h t = 0.070 
for < z < 0.5, Az phot = 0.061 for 0.5 < z < 1, and Az pho , = 0.083 
for 1 < z < 2.5); and (3) the stellar mass of these objects should 
be within 0.1 < M sat /M ccn(ra i < 1.0 for the galaxies in the range 
< z < 2, and within 0.01 < M sat /M ccntral < 1.0 for the galaxies 
in the range < z < 1. An example of satellite galaxy satisfying 
the above criteria is shown in Fig. [2] Finally, we consider different 
redshift bins (see Table [TJ to explore the evolution of the fraction, 
F s;lt , of massive galaxies with satellites. The width of these bins 
were chosen to include a similar number of massive galaxies in 
each bin and have a similar statistics among them. 

We adopted a search radius of 100 kpc. This radius is a com- 
promise between having a large area for finding a significant num- 
ber of satellite candidates that are gravitationally bound to our cen- 
tral massive galaxies but not as large as to be severely contaminated 
by background objects. In any case, we have explored what is the 
effect on our measurements if we select larger radii of exploration. 
We computed the fraction of massive galaxies with satellites for dif- 
ferent search radii CK SC aich=100, 150, 200 and 250 kpc). The results 
of this experiment in the mass range 0.1 < M s;lt /M ccntla i < 1 are 



shown in Fig[3] The numbers presented here are corrected of back- 
ground contamination as it will be explained later. As it is expected, 
we detect an increasing number of massive galaxies with satellites 
as we expand the search radius /^search- The only exception is the 
redshift range 1.1 < z < 2.0 where the fraction of massive galaxies 
with satellites is constant within the error bars. It is worth noting 
that, in general, beyond -K scilr ch = 150 kpc there is not a net increase 
in the fraction of massive galaxies with satellites. Moreover, our re- 
sults are basically unchanged if we use a search radius of 100 kpc 
or 150 kpc. For this reason, on what follows we will present the 
results based on a search radius of 100 kpc as our simulations show 
that this case is affected by the background contamination a factor 
of ~ 2 less than the 150 kpc case. 

Similarly to the selection of the search radius, we have re- 
stricted our potential satellite galaxies to have a redshift difference 
with the central galaxy not larger than the ltr uncertainty in the 
estimation of the photometric redshifts of the Rainbow database. 
Larger redshift differences could be used to include more poten- 
tial candidates but this is transformed also into a larger background 
contribution to our measurements. For instance, we estimated how 
the fraction of massive galaxies with satellites changed when using 
2cr uncertainty in the estimation of the photometric redshifts in- 
stead of lcr. As expected, we found a slight increase (<30%) on the 
fraction but also our error bars increased (by a 50%) by the larger 
amount of background contamination. As these change do not alter 
our main results but increase our error bars, we have used the lcr 
criteria explained in this paper. 

3.1 Background estimation 

Despite we have used photometric redshift information to select our 
potential satellite galaxies, there is still a fraction of objects that 
satisfy the above criteria but are not gravitationally bound to our 
massive galaxies. These objects are counted as satellites because 
the uncertainties on their redshift estimates include them within our 
searching redshift range. These foreground and background objects 
(hereafter we will use the term background to refer to both of them) 
constitute the main source of uncertainty in this kind of studies. 
Consequently, it is key to estimate accurately the background con- 
tamination in order to statistically subtract its contribution from the 
fraction of galaxies with observed satellites. 
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Figure 3. Fraction of massive galaxies having satellites within different 
projected radial distances (search radius, Search) in the mass range 0. 1 < 
^sat/Mcentral < 1.0 for the different redshift bins studied in this work. 



To estimate the fraction of background sources that contami- 
nates our satellite samples we have run a set of simulations. This 
method consists on placing a number of mock massive galaxies 
(equal to the number of our central galaxies) randomly through the 
volume of the catalog. To match our observed redshift distribution 
we assure that in our simulations, the number of mock galaxies 
that are within each redshift bin is the same than in our observed 
sample. Once we have placed our mock galaxies through the cat- 
alog, we count which fraction of these mock galaxies have satel- 
lites around them taking into account our criteria of redshift and 
distances explained above. This procedure is repeated two million 
times to have a robust estimation of the fraction of mock galaxies 
with satellites. We call this average fraction 5 s i mu i. Additionally, 
these simulations allowed us to estimate the scatter in the fraction 
of galaxies that have contaminants. We use this scatter as an esti- 
mation of the error of our real measurements. We consider then this 
fraction to be representative of the background affecting our real 
central sample. The galaxies in our mock samples keep fixed the 
parameters of the massive galaxies (i.e., stellar masses and Sersic 
indices). 

Taking into account that the observed fraction of galaxies with 
satellites, _F obs , is the sum of the fraction of galaxies with real satel- 
lites, F sat , plus the fraction of galaxies which have not satellites 
but are affected by contaminants (1 - F sa i)x5 s i mu i we arrive to the 
following expression: 

^obs ^ simul , , . 

f«* = -rjr • (i) 

The results of our simulations are shown in Table Q] From 
these simulations we see that the fraction S s i m ui of massive galaxies 
we expect to be contaminated by false satellites (using our search- 
ing criteria) is -10% for 0.1 < M S a t /M ccnlra i < 1.0 and -25% for 
0.01 < AWAWa, < 1.0. 



3.1.1 Clustering effects 

It is well known that massive galaxies, particularly in the nearby 
Universe, tend to populate regions which are overdense compared 
to the average density of the Universe. This implies that there is an 
excess of probability (which we will term as clustering) of finding 
galaxies that could be misidentified as satellites of our main tar- 
gets. It is worth noting that this probability excess is not related to 
the accuracy of our redshift estimations. Even with all the redshifts 
measured spectroscopically, the effect of clustering will be equally 



relevant in our estimates as this effect is inherent to our inability 
of measuring real distances but distances inferred by recessional 
velocities. In massive cluster of galaxies, with velocity dispersion 
of - 1000 km s -1 , this will limit our accuracy on estimating real 
galaxy associations. 

Being the clustering a local effect, ideally one would like to 
measure its influence as closer as possible to the central galaxy. 
In practice, this is done by measuring the amount of satellite candi- 
dates in differen t annuli beyond our search radius JChen et alJioOrl 
iLiu et a i] |201 ]]). We will call to the fraction of massive galaxies 
having satellites in these annuli as S cluster- This fraction measures 
both the effect of the background contamination plus the excess 
over this background due to the clustering. This method has the 
disadvantage, compared to the simulations that we have conducted 
above, that is statistically more uncertain. 5 c i„ stcl can be measured 
only around our massive galaxies and this number is relatively 
small. For this reason, S duster is determined with an error larger than 
S simui- We count the satellites in 9 different annuli in the radial range 
100 < R < 330 kpc (the size of each annuli was selected to contain 
the same area than the searching area within 100 kpc). We find, as 
expected, that the number of satellites decreases in the outer an- 
nuli, reaching asymptotically (within the errors) the values we get 
using the first background estimation method. However, in general, 
and particularly for the lower redshift bins, the number of detected 
satellites is higher for the inner annuli than in the random case, and 
therefore, the clustering is not negligible. As we noted before, the 
detection of satellites does not increase at R xarc h > 150 kpc. For 
this reason, and as a compromise between proximity to the mas- 
sive galaxies and having enough statistics, we have used the aver- 
age detections of satellites in the two annuli closer to R=150 kpc 
(173 < R < 200 kpc and 200 < R < 224 kpc) to estimate the 
effect of the clustering. The uncertainty at measuring 5 cluster is not 
straightforward to calculate and we have decided to estimate that 
value summing quadratically the background uncertainty measured 
in the simulations estimating 5 s imui plus the dispersion between the 
two different radial annuli used in the clustering determination. 

The significance of the clustering is quantified in Table [T] We 
find that above z > 1 the clustering is playing a minor role as Setter 
and 5" s i mll i are very much alike within the errors. However, at z < 1, 
^cluster - 1.5 xSsimut. As it is expected, the effect of the clustering is 
more relevant at lower redshifts. At high redshifts, the overdensities 
are less significant as the large scale structures would be still not 
completely formed. 



4 RESULTS 

In Table Q] we summarize the results obtained in this work. For 
each redshift bin, we present the fraction of galaxies with satel- 
lites initially found for our sample of massive galaxies, F D b s , the 
background estimate S s j mu i derived from the mock catalogues, and 
the final fraction of massive galaxies with satellites, F sal , after the 
correction of the background contamination with Eq.Q] The associ- 
ated errors correspond to the standard deviation from the measure- 
ments obtained in the mock catalogues as explained in the previous 
section. In addition, we include the expected contamination due to 
the clustering estimate, S duster, and the fraction of massive galaxies 
with satellites after this correction F c i ustor . 

Our results are also illustrated in Fig. [4] Our main result is seen 
in the upper left panel of Fig. [4] the fraction of massive galaxies 
with satellites, within a projected radial distance of 100 kpc, in the 
range 0. 1 < M S a t /M ccnll ai < 1 remains basically constant (17 ± 3%) 
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Figure 4. Fraction of massive galaxies having satellites (and their properties) within a projected radial distance of 100 kpc for different redshift bins. Upper 
panels show the fraction of massive galaxies with satellites in the mass range 0. 1 < M sat /M centra i < 1.0 (orange dots) and the fraction of massive galaxies 
with satellites in the mass range 0.01 < M sat /M ccntra i < 1.0 (magenta dots). In the bottom panels, we explore the fraction of massive galaxies with satellites 
in the mass range 0. 1 < M sat /M ccntra i < 1.0 when our sample is split depending on the Sersic index (morphology) of the galaxies (red dots for n > 2.5 
-spheroids- and blue stars for n < 2.5 -disks). The horizontal dotted lines in the upper left panel correspond to the average fitted values to our findings. Central 
panels: Mean projected distances where the satellites are found after statistical correction of the background. The dashed line indicates the average projected 
distance obtained for mock satellites in the simulations (i.e. this is the expected distance if the satellites where an artifact of the background contamination). 
Right panels: Mean mass ratios between the satellites and their massive galaxies. The horizontal bars indicate the range of redshifts considered for each 
measurement. For clarity, we have slightly shifted the data corresponding to the spheroids-like objects. 



in the redshift interval < z < 2. To have a z = comparison, we 
have added the measurement from lLiu et all ( 1201 1) using the SDSS 
sample. They find that at z = the fraction of massive galaxies 
with satellites in the mass range and projected radius explored here 
is very similar. In the same panel, we show the same analysis up 
to z = 1 for satellite galaxies with 0.01 < M sat /M ccntra | < 1.0. 
Although a little bit noisier due to the lower statistics, our findings 
agree with a relative constant fraction (3 1 ±6%) of massive galaxies 
having such type of satellites. 

Our sample of massive galaxies is large enough that we can 
explore whether the fraction of massive galaxies with satellites in 
the range 0.1 < M sat /M ccntra i < 1 depends on the morphology of 
the massive galaxy. We have used the Sersic index as a proxy to the 
galaxy morphology. In the nearby universe, galaxies with n < 2.5 
are mostly disc-like o bjects, whereas galaxies with n > 2.5 are 
mainly spheroids (e.g., Andredakis et al.ll995l ; lBlanton et alj2003l ; 
iRavindranath et al.ll2004h . We have used the published Sersic in- 
dexes provided by T07 to separate our galaxies. We illustrate our 
results in the bottom left panel of Fig. [4] There is a hint that massive 
galaxies with spheroid-like morphologies tend to have a larger frac- 
tion (a factor of 2-3) of galaxies with satellites than disk-like mas- 
sive objects. This result is more prominent at low redshift where the 
clustering of the massive spheroid population could be an issue. 

We can repeat the same exercise but this time using the clus- 
tering correction (which contains also the background effect) to ex- 
plore how our results depend on this effect. The comparison be- 
tween the two types of corrections are shown in Fig. [5] In general, 



the correction due to the clustering decreases the fraction of mas- 
sive galaxies which contains satellites. This is now 12 ± 2% in the 
redshift interval < z < 2 for 0.1 < M sat /M centra i < 1 and 23 + 4% 
for 0.01 < M sat /M cenlra , < 1.0 up to z = 1. 

4.1 Robustness of the results 

The results presented in this paper are the product of combining two 
different datasets: the T07 sample of massive galaxies and the Rain- 
bow catalogue. In addition, we have used photometric redshifts 
and, when available (54% of the time), spectroscopic redshifts. We 
have checked how robust are our results to the use of more ho- 
mogeneous dataset, as using only a sample of massive galaxies 
with spectroscopic redshifts or, to base our full analysis only on 
the Rainbow database. 

In our first test, we used only the central galaxies in our sample 
with spectroscopic redshifts and counted which fraction has satel- 
lites following the same procedure explained above and after cor- 
recting for the background. The width of the redshift bins for this 
sample were again chosen to include a similar number of massive 
galaxies and have a similar statistics to compare with the other sam- 
ples. The output of this test is shown in Fig(6] We get, for the case of 
0.1 < M sat /M ccntla | < 1 up to z = 2, an average fraction of 19 ± 4%. 
It can be seen that this result is in full agreement with our previ- 
ous estimation for this quantity. In a second test, we have taken 
the redshifts and stellar masses only from the Rainbow catalog to 
check whether there are systematic effects due to the use of com- 
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Table 1. Fraction of massive galaxies with satellites at different redshifts. For each redshift range we present the number of massive galaxies A^nu-al m eacn 
bin (number of galaxies with spectroscopic redshifts in brackets), the observed fraction of massive galaxies with satellites F b s , the estimate of the background 
contamination S s j ml ,i, and the estimate of the clustering effect S cluster- Finally, we present the final fraction of massive galaxies with satellites when i) the 
correction for the background contamination (f sat ) or ii) the clustering effect (F c i uster ) is applied. 



Redshift range 


"central 


Fobs 


§ simul 


S cluster 


fsat 


^cluster 




(N with spec z) 












All galaxies 


0.10 < M sa t/M ccntral < 1.00 












0.20 < z < 0.75 


197 (130) 


0.29 


0.09 ± 0.02 


0.15 ±0.02 


0.22 ± 0.02 


0.17 ±0.02 


0.75 < z < 0.90 


129 (76) 


0.24 


0.10 ±0.03 


0.17 ±0.03 


0.16 ±0.03 


0.08 ± 0.03 


0.90 < z < 1.10 


142 (99) 


0.25 


0.08 ± 0.02 


0.13 ±0.03 


0.18 ±0.02 


0.13 ±0.02 


1.10 <z< 2.00 


161 (55) 


0.18 


0.09 ± 0.02 


0.09 ± 0.03 


0.10 ±0.02 


0.10 ±0.02 



Q.01 < Msat/AWa, < 1.00 



0.20 <z< 0.55 51 (40) 
0.55 < z < 0.73 70 (42) 
0.73 <£< 1.10 73 (53) 


0.37 
0.53 
0.52 


0.20 ± 0.06 
0.24 ± 0.05 
0.27 ± 0.05 


0.24 ± 0.06 
0.36 ± 0.05 
0.36 ±0.06 


0.22 ± 0.06 
0.38 ± 0.05 
0.34 ± 0.05 


0.16 ±0.06 
0.28 ± 0.05 
0.25 ± 0.06 


Spheroid like (n > 2.5) galaxies 


0.10 < Msat/AWal < 100 
0.20 < z < 0.75 137 
0.75 <z< 1.10 176 
1.10 <z< 2.00 85 


0.34 
0.27 
0.18 


0.09 ± 0.03 
0.08 ± 0.02 
0.08 ± 0.03 


0.16 ±0.03 
0.14 ±0.04 
0.08 ± 0.03 


0.28 ± 0.03 
0.21 ±0.02 
0.10 ±0.03 


0.23 ± 0.03 
0.15 ±0.04 
0.10 ±0.03 


Disk like (n < 2.5) galaxies 


0.10<M sat /M centn i < 1.00 
0.20 < z < 0.75 60 
0.75 <z< 1.10 95 
1.10 <z< 2.00 76 


0.18 
0.19 
0.18 


0.10 ±0.03 
0.09 ± 0.03 
0.09 ± 0.04 


0.15 ±0.04 
0.16 ±0.04 
0.09 ± 0.05 


0.09 ± 0.03 
0.10 ±0.04 
0.10 ±0.04 


0.05 ± 0.04 
0.04 ± 0.04 
0.11 ±0.05 



billing different samples. We get, in this case, a fraction of 17 ± 2%. 
Again, this result agrees perfectly with the original estimate. We 
conclude, accordingly, that our results are robust to both the use of 
spectroscopic redshifts only and to the mixing of different datasets. 

Another test that we have conducted is to check whether our 
results are robust to a change in the stellar mass limit at selecting 
our massive galaxies. As it is illustrated in Fig. [T] the galaxies at 
higher redshifts are slightly more massive than the bulk of objects 
at lower redshift to guarantee that we can study satellites above a 
given mass ratio along the full redshift range. We have checked how 
our results change if we select only massive galaxies with stellar 
masses above 2x10" M G . We have done this exercise for the case 
0.1 < M sat /M ccntral < 1 up to z = 2. With the new mass limit, 
we get an average fraction of ~ 23%. This is in good agreement 
with our original estimation for this fraction. Again, increasing the 
stellar mass limit does not alter substantially the fraction of massive 
galaxies with redshift. 



4.2 Properties of the satellite galaxies 

In addition to counting which fraction of the massive galaxies have 
satellites, we can also estimate the average projected radial dis- 
tances of these satellites and the average mass ratios between the 
satellites and the massive objects. To estimate these quantities prop- 
erly, we need to correct statistically by the effect of the contami- 
nants. This can be done using the following expression: 

, , *ob>i S simul , , 

< gsat >= -T— < gobs > ~ < gsimul > (2) 

^sat ^sat 



where < <2 b s > is the observed mean value of the property Q (i.e. 
the projected radial distance or the mass ratio), < <2 s j mm > is the 
mean value obtained from the mock massive galaxies (i.e. the val- 
ues that are found for the contaminants) and < Q SM > is the value 
after the correction. 

The mean projected radial projected distances (in kpc) of the 
satellites and their mean mass ratios M s;lt /M cen t r ai are compiled in 
Tableland shown in Fig. [4] We plot with a dashed line the aver- 
age projected distance (~ 72 kpc) of the background galaxies de- 
tected as fake satellites in the simulations. After correcting by the 
effect of the contaminants, we find that our satellite galaxies are at a 
typical projected radial distance of ~40 kpc. This value is well be- 
low the expectation from a random distribution, suggesting that the 
satellites are gravitationally bounded to their central galaxies. This 
average distance seems to be pretty much independent (within the 
errors) of the satellite mass, the morphological type of the central 
galaxy and the redshift of the system. 

Finally, we show in the right panels of Fig. f4]the mean mass 
ratio M sat /M ccntral in each redshift bin after the statistical correction. 
We find that M sa( is around 0.36 M centla i when we explore satellite 
galaxies within a mass ratio of 0.1 < M sat /M ct . utra ] < 1 (this value is 
0.28 when we use the clustering correction). If we explore down to 
a mass ratio of 0.01 then M sM is around 0.15 M ccnlia i (0.14 correct- 
ing by the clustering effect). It is worth noting that in both cases, the 
mean masses of our satellites are over 10 10 M o and consequently we 
are detecting satellites with large masses. When we split the sam- 
ple depending on their Sersic indices (bottom right panel in Fig.[4]l, 
there are not significant differences, within the errors, between both 
samples. Again we find that the satellites of massive galaxies are 
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1 : 1 0— background O 
1 : 1 0— clustering • 
1:100-background O 




I i i i i I i i i i I i i i i I 

0.5 1 1.5 



Table 2. Mean projected radial distances between the satellites and the cen- 
tral galaxies and their mean mass ratios (M sat /M ccntra i). These values corre- 
spond to the case where only the background correction has been applied. 



Redshift range 


Radial Distance 


AWMcentral 




(kpc) 




All galaxies 


0.10 < M sat /M ccnCral < 1.00 




0.20 < z < 0.75 


41 ±4 


0.35 ± 0.03 


0.75 <z< 0.90 


38 ±6 


0.40 ± 0.07 


0.90 < z < 1.10 


49 ±6 


0.36 ± 0.05 


1.10 < z < 2.00 


28 ±6 


0.34 ± 0.08 



0.01 < M sat /M t , nCral < 1.00 

0.20 < z < 0.55 25 ± 7 0.23 ± 0.07 

0.55 <z< 0.73 41 ±5 0.16 ±0.02 

0.73 <z< 1.00 23 ±3 0.19 ±0.03 



Spheroid - like (n > 2.5) galaxies 

0.10 <M sat /M ccntral < 1.00 

0.20 < z < 0.75 42 ± 4 0.33 ± 0.03 

0.75 <z< 1.10 47 ±3 0.38 ±0.05 

1.10 <z< 2.00 42 ±11 0.25 ±0.06 



Disk - like (n < 2.5) galaxies 



a CN 
b_ O 



0.5 



spheroids — background 
spheroids— clustering 
disks-background 
disks — clustering 



1.5 



Figure 5. Fraction of massive galaxies with satellites after correction with 
the background contribution computed from the simulations (open sym- 
bols) and when the clustering estimation is used (filled symbols). Top panel: 
The fraction of massive galaxies with satellites is separated in two groups: 
0.1 < M sat /M ceIlt ral < 1.0 (blue triangles) and 0.01 < M sat /M centra i < 1.0 
(red squares). Bottom panel: Fraction of massive galaxies in the mass range 
0.1 < A/ sal /M ci:Iltra i < 1.0 depending on the Sersic index (morphology) of 
the central galaxies. 



0.10 < M sat /M ccnlral < 1.00 

0.20 <z< 0.75 33 ± 12 0.33 ±0.13 

0.75 <z< 1.10 33 ± 13 0.30 ±0.11 

1.10 <z< 2.00 19 ± 7 0.44 ±0.16 



similar in their mean properties (distances and mass ratios) inde- 
pendently of the morphological type of their central galaxies and 
redshift. 



1:10 mixed • 
1:10 spec O 
:10 Rainbow d 



0.5 



1.5 



Figure 6. Testing the robustness of the fraction of massive galaxies with 
satellites in the mass range 0. 1 < A/ sat /iW c( . ntla i < 1.0 for different datasets: 
our original sample described in Sect. 2 (red points), a purely spectroscopic 
sample (magenta circles) and a sample based only on the Rainbow cata- 
logue (blue squares). 



5 DISCUSSION AND CONCLUSIONS 

The results of this paper support a picture where the fraction of 
massive (M star ~ 10" M Q ) galaxies with satellites, within a pro- 
jected radius of 100 kpc, has not changed with time since z ~ 2. 
This fraction remains around ~ 15 % for galaxies with satellites 
with mass M star > 10 10 M G and around ~ 30% if we explore satel- 
lites with masses M star > 10 9 M G up to z = 1. 

Interestingly, we find a hint that the fraction of massive galax- 
ies with satellites is larger (a factor of around 2 to 3) for those galax- 
ies with spheroid-like morphologies than for galaxies with disk-like 
appearance (much evident for z < 1.1). This fact could be linked to 
the different size growth we observe for these two types of objects 
as cosmic time increases. In fact, spheroid galaxies are known for 
growing more dramatically in size since z ~ 3 than disk galaxies 
(see e.g., T07 iBuitrago et aUl2008h . It could be also possible that 
the fraction of satellites difference between spheroid and disk-like 
galaxies would be just an effect of the clustering, more relevant at 
lower redshifts (Sect. [3~TTTV However, this difference remains even 
when this effect is taken into account (Fig. |5). We remark, how- 
ever, that it is difficult to correct the clustering effect accurately. 
With the present dataset, a mild redshift evolution of the fraction of 
spheroid-like galaxies with satellites can not be excluded. 

Due to the enorm ous uncertainty on the merging timescales 
(e.g., lLotz et al]|201ll) . it is beyond the scope of this work to 
estimate a robust merger rate associated to our measurements. 
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Nonetheless, we can make a crude estimation on the number of 
mergers a massive galaxy experiences since a given z according 
to the following expression: N m = T(z) x F sal /r m , where T(z) 
is the interval of cosmic time since a given z to now, and r m is 
the merging timescale of the satellite within a given radius. For 
each massive g alaxy at z = 2 and assuming T m ~1,5 Gyr (e.g. 
lLotz et al .2011), we would expect that the number of mergers with 
mass ratio around 1:3 would be ~1 (~2 in the case of an spheroid- 
like galaxy) since that epoch. For a massive galaxy at z = 1, we 
would expect that the number of mergers with mass ratio around 
1:6 would be ~ 1.5 since that redshift. Again, these numbers are un- 
certain and very much dependent on the exact merging timescale 
which is a function of the baryo nic mass ratio and the model used 
to estimate this quant i ty (e.g.. iBluck et al J [20091 : IConselice et al.l 
120091 : lLotz et alj|201 ll : flVIan et alj|201 ll) . These numbers of merg- 
ers, however, are slightly lower (although the exact amount is dif- 
ficult to quantify) than the expected number of mergers obtained 
using theor etical recipes for th e size incre ase of a galaxy a fter a 
merger (see iTruiillo et al .11201 ll) . Recently, lijiuck et al.l jioTll) find 
that a massive galaxy (M star > 10 n M G ) will experience on aver- 
age N,„ = (1.1 ± 0.2)/r m minor mergers over the redshift range 
z = 1.7-3. This would mean a final A',,, ~ 1 using the r m consid- 
ered in our work. If this is confirmed, it will point out to the pos- 
sibility that the merging activity at those redshifts would be higher 
than at lower redshifts. When extrapolated in redshift, they find a 
total final number of minor mergers of A',,, = (4.5 ± 2.9)/r m from 
z = 3, although once more, the large errors make very uncertain to 
constrain the final number of experienced mergers. 

At present, there are a few cosmological simulations where the 
size growth of the massive galaxies is explained by the accretion of 
minor satellites (see lNaab et al.l2009l : rOser et al.l201 ll) . It would be 
straightforward to compare our findings with those cosmological 
simulations of galaxy formation and check whether the fractions 
that we find are recovered in such theoretical analysis. If this was 
the case, the support to the minor merging mechanism as the main 
responsible for the size evolution of the massive galaxies will be 
greatly enhanced. 
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